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A simplified thermodynamic approach of the incompressible axisymmetric Euler equations is con- 
sidered based on the conservation of helicity, angular momentum and microscopic energy. Statistical 
equilibrium states are obtained by maximizing the Boltzmann entropy under these sole constraints. 
We assume that these constraints are selected by the properties of forcing and dissipation. The fluc- 
tuations are found to be Gaussian while the mean flow is in a Beltrami state. Furthermore, we show 
that the maximization of entropy at fixed helicity, angular momentum and microscopic energy is 
equivalent to the minimization of macroscopic energy at fixed helicity and angular momentum. This 
provides a justification of this selective decay principle from statistical mechanics. These theoretical 
predictions are in good agreement with experiments of a von Karman turbulent flow and provide a 
way to measure the temperature of turbulence and check Fluctuation-Dissipation Relations (FDR) . 
Relaxation equations are derived that could provide an effective description of the dynamics towards 
the Beltrami state and the progressive emergence of a Gaussian distribution. They can also provide 
a numerical algorithm to determine maximum entropy states or minimum energy states. 



PACS numbers: 



I. INTRODUCTION 



In a turbulent flow, the number of degrees of freedom 
scales like i?e^/^, where Re is the Reynolds number, and 
can reach 10^^ for atmospheric-like flows, comparable to 
the Avogadro number. This is beyond the present ca- 
pacity of computers. For example, numerical simulations 
of a von Karman (VK) turbulent flow at Re = 10^, a 
standard laboratory flow used for turbulence studies (see 
below), would require resolutions of the order of 10^^ grid 
points and integration times of the order 10^ years of cpu 
with current computers. This conclusion justifies the in- 
troduction of turbulence models to reduce the number 
of degrees of freedom and make turbulence amenable to 
numerical simulation or theoretical understanding. This 
goal cannot be reached unless the different components 
of turbulence and their interactions are identified. 

Turbulence being intrinsically a stochastic process, it 
can be decomposed in two components: the mean flow 
and the fluctuations around it. A good turbulence model 
should therefore be able to predict both the structure of 
the mean flow, and its influence on and through fluctu- 
ations, within a reduced number of degrees of freedom. 
This kind of information is typically provided by statisti- 
cal mechanics. Can we adapt statistical methods to deal 
with the turbulence problem? 
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This program has been pioneered by Onsager 
Montgomery & Joyce and Lundgren & Pointin Isi in 
the framework of two-dimensional point vortices. In the 
last decade, this statistical approach has been extended 
by Miller Q and Robert and Sommeria Q to simplified 
2D or quasi 2D flows with continuous vorticity. Even 
more recently, Leprovost et al. Q have shown that the 
2D formalism could actually be applied to a typical 2D 
1/2 situation, an axisymmetric flow. They obtained a re- 
lationship that gives the general shape of stationary solu- 
tions (mean flows) of the axisymmetric Euler equations. 
This relationship has been tested and confirmed experi- 
mentally in a turbulent von Karman flow by Monchaux 
et al. Q who observed that, at high Reynolds numbers, 
the selected shape is Beltrami, with vorticity and velocity 
aligned everywhere. As discussed in Appendix I A 11 such 
a shape cannot be obtained with the thermodynamical 
approach of Leprovost et al. Q • In the present work, we 
revisit the theoretical tools in order to capture Beltrami 
states as statistical equilibrium states. We also extend 
the computations one step further by considering fluctu- 
ations around mean field. 

Specifically, we develop a simplified thermodynamic 
approach based on the conservation of helicity, angu- 
lar momentum and microscopic energy. We assume that 
these constraints are selected by the properties of forcing 
and dissipation. From a maximum entropy principle we 
derive the mean flow and the fluctuations around it. We 
find that the mean flow is in a Beltrami state and that the 
fluctuations are Gaussian. We also show that the maxi- 
mization of entropy at fixed helicity, angular momentum 
and microscopic energy is equivalent to the minimization 
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of macroscopic energy at fixed helicity and angular mo- 
mentum. This justifies from statistical mechanics a se- 
lective decay principle introduced previously from phe- 
nomenological arguments @ . We use the mean field the- 
ory to link the fluctuations to the response of the mean 
flow to perturbations (susceptibility) and to the temper- 
ature in a way reminiscent of the Fluctuation-Dissipation 
Theorem. This provides a way to measure the temper- 
ature of turbulence through the fluctuation level. The 
analogy with 2D turbulence is discussed. In fact, due to 
the dual nature of axisymmetric flows, intermediate be- 
tween 2D and 3D turbulence, we find the emergence of 
two different effective temperatures in the fluctuations. 
One temperature, characterizing velocity fluctuations, is 
related to the formation of coherent structures like in 2D 
turbulence. Another one, characterizing vorticity fluctu- 
ations, is related to 3D vorticity stretching and diverges 
with increasing resolution. These predictions have been 
tested in companion papers [l^ based on PIV mea- 
surements in a turbulent von Karman flow and are in 
fair agreement with observations. 

The paper is organized as follows. In Section |lll wc 
recall the basic problematics associated with statistical 
mechanics of turbulence, and formulate our hypotheses 
and the associated theoretical framework. In Section Hill 
we recall the stationary solutions and the conservation 
laws (the backbone of the statistical mechanics approach) 
of axisymmetric flows. In Sec. IIV[ we recall the phe- 
nomenological selective decay principle leading to Bel- 
trami flows. Section |V] is devoted to the computation of 
the statistical equilibrium states of axisymmetric inviscid 
flows using mean field theory. We derive the Gibbs states 
and the fluctuation-dissipation relations (FDR) with two 
different mean field approximations. In each case, the 
mean flow is in a Beltrami state and the fluctuations 
are Gaussian. In Section I VII we make the connection 
between different variational principles that character- 
ize the equilibrium states. For each principle, we pro- 
pose a set of relaxation equations that can be used as 
a numerical algorithm to solve the variational problem. 
These relaxation equations can also provide an effective 
description of the relaxation of the system towards maxi- 
mum entropy states. Finally, we justify through statisti- 
cal mechanics the phenomcnological principle according 
to which: "the mean flow should minimize the macro- 
scopic energy at fixed helicity and angular momentum" . 

II. HYPOTHESES AND THEORETICAL 
FRAMEWORK 

A. Turbulence, Navier-Stokes equations and 
classical statistical mechanics 

A turbulent flow is described by the Navier-Stokes 
equations 

^ + u-Vu = Vp + i/V^u + Z, (1) 

ot p 



where u is the velocity, p the pressure, p the fluid density, 
v its kinematic viscosity, and / a forcing. In the absence 
of forcing, the velocity decays to zero due to the dissipa- 
tion, so that turbulence is an intrinsic out-of-equilibrium 
problem. In the sequel, we focus on the simplest situa- 
tion, where forcing and dissipation equilibrate on aver- 
age, so that stationary states can arise. The goal of the 
present paper is to describe these stationary states and 
the fluctuations around them using tools borrowed from 
classical statistical mechanics. Speciflcally, we are going 
to introduce a Hamiltonian system, perform equilibrium 
or near equilibrium statistical mechanics and compute its 
equilibrium states. 



B. Stationary Navier-Stokes solutions vs solutions 
of Euler equations 

Since forcing and dissipation equilibrate on average 
for stationary solutions of the Navier-Stokes equations, 
it seems natural to consider this limiting case flrst in our 
quest of a framework suitable for classical statistical me- 
chanics. In such a limit, we get the Euler equations 

-^+uVu = — Vp. (2) 

ot p 

This is indeed a Hamiltonian system as long as one con- 
siders regular solutions such as those based on finite 
Galcrkin expansions. In 2D turbulence, the considera- 
tion of Euler solutions to describe Navier-Stokes station- 
ary solutions is well accepted, based on the remark that 
the vorticity cannot blow up and that the limit v 
is usually well behaved under reasonable regularity hy- 
pothesis. In 3D turbulence, this hypothesis is still con- 
troversial since Onsager [l[. One major problem is that 
one cannot exclude vorticity blows up in 3D, that would 
make the limit — >■ singular. A signature of this ef- 
fect is the famous 4/5 law of homogeneous turbulence 
that links energy dissipation to the third moment of the 
velocity increments, independent of any viscosity. For 
this reason, Onsager [l| suggested to consider weak solu- 
tions of the Euler equations to describe stationary states 
of the Navier-Stokes equations, thereby allowing a finite 
amount of energy dissipation even in the absence of vis- 
cosity. This suggestion was developed recently in an el- 
egant way by Duchon and Robert [ll[. However, weak 
solutions are not directly amenable to methods of classi- 
cal statistical physics, and nobody has yet succeeded to 
follow to the end Onsager's suggestion. In the present 
case, we overcome this difficulty by considering only reg- 
ular solutions of the Euler equation. The bonus is that 
we deal with a Hamiltonian system to which we can ap- 
ply statistical mechanics. The malus is that we may have 
lost any connection with actual turbulence. However, in 
companion papers 0, |^, we compare our theoretical 
predictions with actual experimental turbulent flows and 
show that they basically agree. In other words, it seems 
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that stationary states and fluctuations of an out-of equi- 
librium system, the forced Navier-Stokes equations, can 
be described by statistical equilibrium states and fluctua- 
tions of the Euler equations without forcing and dissipa- 
tion. 



C. The Euler system and conservation laws 

The Euler equations for regular solutions are charac- 
terized by a number of conservation laws that depend on 
the geometry and on the dimension of the system. In 2D 
turbulence, for example, the conservation laws arc the 
kinetic energy E ~ ^ J u^dx, the enstrophy ^ = J uj'^dx, 
where wz = V x u is the vorticity, and, more generally, 
any function of the vorticity (Casimirs). In 3D turbu- 
lence, the generic conserved quantities are the kinetic 
energy E and the helicity H = J u ■ oj dx. Additional 
conservation laws are possible in the presence of addi- 
tional symmetries, such as axisymmetry, see 0]. 

In the presence of forcing and dissipation, these conser- 
vation laws are altered. In the sequel, we shall postulate 
that the balance between forcing and dissipation selects 
some particularly relevant conservation laws among the 
infinity of inviscid invariants. In particular, we shall ar- 
gue that there exists relevant situations in which the only 
conserved quantities are the microscopic energy E and 
the helicity H. This property holds, for example, for very 
simple solutions of the Euler equations such that the ve- 
locity and the vorticity are aligned everywhere in the flow 
(Beltrami state). Our aim in this paper is not to deter- 
mine the mechanisms that select these invariants. This 
is a complicated problem that depends on the properties 
of forcing and dissipation and on the Reynolds number. 
However, to motivate our approach, we show in compan- 
ion papers 0, |^, that our assumptions are consistent 
with experimental results in the limit of large Reynolds 
numbers. 



D. Boundary conditions 

When we try to model a turbulent flow by numerical 
simulations, the choice of boundary conditions is a ma- 
jor and complicated problem. Indeed, the results turn 
out to depend sensitively upon what boundary conditions 
have been chosen. For example, for two-dimensional 
Navier-Stokes turbulence, qualitatively different results 
have been obtained for free-slip boundary conditions in 
rigid squares and rectangles 0], for rectangu- 

lar periodic boundaries [isl [l6| . for no-slip circular 
boundaries [13, [I!] , for free-slip circular boundaries [l^ , 
for stress- free circular boundaries [ll], and for no-slip 
square boundaries For our problem, we argue that 
numerical simulations cannot deal with sufficiently high 
Reynolds numbers (see the Introduction) and we rather 
focus on an experimental device, namely a von Karman 
flow 0, H- In. that case, the boundary conditions are 



automatically determined by the experimental geometry 
and the forcing. 

Since our theoretical approach is inviscid, we cannot 
hope to describe the experimental flow close to its bound- 
aries. We rather claim that, when forcing and dissipation 
equilibrate each other, this flow can be reasonably de- 
scribed with an inviscid approach far enough from the 
boundaries. Our theory will be derived with the inviscid 
boundary conditions reflecting a cylindrical closed do- 
main (i.e. tp = on the boundary). The comparison 
with experiments, detailed in [lo| . will be performed over 
a sub-domain of the experiment, far from boundary and 
forcing. As already shown in 0, @1 j this restriction is suf- 
ficient to obtain good agreement between inviscid theory 
and experimental fields, at the zeroth order approxima- 
tion (i.e. if one considers the mean fiow topology and its 
small fluctuations). 



III. EULER EQUATION IN THE 
AXISYMMETRIC CASE 



A convenient formulation of the axisymmetric 
Euler equations 



In the axisymmetric case, the incompressible Euler 
equations take the form 
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1 dp 
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(3) 

(4) 

(5) 
(6) 



where (m^, ug,Uz) are the velocity components in a cylin- 
drical referential (r, 0, z). Here, r runs from to i? and z 
from to 2h (we take the origin of the z axis at the bot- 
tom of the domain). Furthermore, we choose the length 
unit such that the total volume is unity: J rdrdz = 1 
(due to the axial symmetry, we systematically divide all 
the volume integrals by 27r). It was shown in that 
the axisymmetric incompressible Euler equations can be 
rewritten in a simplified form in terms of ct, ^ and ip, 
where a = rug is the angular momentum, ^ is the po- 
tential vorticity related to the azimuthal component of 
the vorticity by ^ = (-L>g/r, and tp is the streamfunction 
associated with the poloidal component of the velocity: 



U = Ugeg + V X —eg 



(7) 
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Note that Ur = —dzip/f and Uz = drip/r. The axisym- 
metric Euler equations can then be recast as @ : 



da 
'dt 



+ {V,a} = 0, 



1 d^^/j 
2y dz^ 



(8) 
(9) 

(10) 



where y = { , } is the Poisson bracket {{ijj^ff)} = 

dyipdz4> — dz'ipdy'P) a-nd A* is a pseudo-Laplacian. 

A few general remarks are in order regarding this spe- 
cial case: i) One sees from Eq. that the angular mo- 
mentum is conserved by the fluid particles and can only 
be mixed through the Euler dynamics. This is the ana- 
log of vorticity mixing in 2D turbulence and it justifies 
the introduction of a mixing entropy for the distribution 
of angular momentum (see Sec. IVB[) . We will see that 
we can establish a close parallel with the statistical me- 
chanics of 2D turbulence to determine the distribution of 
angular momentum, ii) By contrast, the potential vor- 
ticity is stirred like in 3D turbulence and not conserved. 
Therefore, the distribution of potential vorticity is more 
difficult to investigate and this can lead to complicated 
problems such as cascade towards small scales, forma- 
tion of singularities etc. In fact, we shall show in the 
companion paper [2l| that the statistical theory predicts 
the existence of large-scale coherent structures (like in 
2D) but that these states are unstable saddle points and 
should cascade towards smaller and smaller scales (like 
in 3D). However, these large-scale structures can have a 
very long lifetime because, being saddle points of entropy, 
they are unstable only for some particular perturbations. 
If the dynamics does not spontaneously generate these 
optimal perturbations, the system can remain "frozen" in 
a saddle point of entropy for a long time [111, . There- 
fore, we are truly in a situation intermediate between 2D 
turbulence and 3D turbulence. 



B. Stationary states 

The axisymmetric Euler equations admit an infinite 
number of steady states. The general form of stationary 
solutions of the axisymmetric Euler equations (|51 fTU|) has 
been established in Q. They are given by 



^ = /(V'), (11) 

C = — ^ — + gW, (12) 



where / and g are arbitrary functions. When / is linear 
/(x) = Xx and g = 0, the vorticity and velocity are 
aligned everywhere a; = Au and the stationary fiow is a 
Beltrami state. An alternative form also established in 



^R'i'T) - J = Q(a), 



(13) 
(14) 



where R and Q are arbitrary functions. The relation be- 
tween Eqs. pllin|) and Eqs. ([TMI)) is developed in ^ 
provided that some invertibility properties for the func- 
tions are assumed. 



C. Conservation laws 

Axisymmetric inviscid flows admit an infinite number 
of conserved quantities, namely the total energy 



E 



1 



{ul + ul + ul)rdrdz 



= ^ y i-ipdydz + ^ y ^dydz, 

the Casimirs 

Jg = y G{a)dydz, 
and the generalized helicities 



Hf = j iF{a)dydz, 



(15) 



(16) 



(17) 



where G and F are any (regular) functions. In the se- 
quel, we also introduce the notation _ff„ and /„ for the 
case where F ov G are power-laws z". In particular, / ~ 
Ii = J adydz (angular momentum), F = Hq = J ^dydz 
(circulation) and H = Hi = J ^adydz (helicity) are con- 
served. 



D. Energy-helicity-Casimir functional 

From the integral constraints discussed previously, a 
generalization of the Arnol'd energy-Casimir functional 
has also been introduced in Q. This is the energy- 
helicity-Casimir functional A = E + Iq + Hp. Consider 
the optimization problem 



min / max{ ct] }. 



(18) 



A critical point of this functional determines a steady 
state of the axisymmetric equations. Indeed, writing 



6A ^ 5{E + Ig + Hp) =0, 
and taking variations on cr and ^, we obtain 

^ + F{a)=0, 



2y 



G'{a)+iF'{a)^0, 



(19) 



(20) 



(21) 
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and we recover the equations (|13m4p characterizing a 
steady solution of the axisymmetric Euler equations. The 
fact that we obtain all the steady states means that the 
quantities given by Eqs. (|15ll6ll7p are the unique invari- 
ants of the axisymmetric incompressible Euler equations 
0. Furthermore, if the critical point of Eq. (|T8l) is a 
maximum or a minimum of A then this steady state is 
nonlinearly dynamically stable. In many cases, we shall 
restrict ourselves to formal nonlinear stability [2^. We 
consider small perturbations and we only require that 
the critical point is a (local) maximum or minimum of A 
such that the second order variations 



S^A = 



1 



S(,Sip dyd 
1 



{5a) 



■ dydz 



Ay 

F"(a)^(5af dydz 
F'{a)S^Sadydz, (22) 



2 G"{a)iSafdydz+ ^ 



are definite positive or definite negative for all perturba- 
tions 6a and 6^. Formal stability implies linear stability 
(in that case 6^ A can be used as a norm) but it does not 
irn|)ly nonlinear stability for infinite dimensional systems 

On the other hand, the optimization problem given by 
Eq. (|18p provides just a sufficient condition of nonlinear 
dynamical stability. More refined stability conditions can 
be obtained by adding some constraints in the optimiza- 
tion problem [2^ For example, the minimization 
problem 

inin{ E[t a] \ Hf[^, a] = Hp, Ig[L ^] - Ig }, (23) 



is more refined than 
min{ ct] - 

?.cr 



(24) 



in the sense that a solution of Eq. p4)) is always a solu- 
tion of the more constrained problem given by Eq. (|23p . 
but the reciprocal may be wrong. This is similar to en- 
semble inequivalence in statistical mechanics where dif- 
ferent ensembles have the same critical points but not 
necessarily the same maxima or minima (giving rise to 
different stability criteria) [1^ [l^ . Ensemble inequiva- 
lence is generic for systems with long-range interactions 
like turbulence. 



The time variations of A are given by 

A = j X {tlj + F{a))dydz 

+ / ^ ( J + G'(a) + (F'{a)^ dydz. (26) 

To determine the functions X and we maximize the 
rate of production (resp. dissipation) of A with the con- 
straints 



Y2 y2 



(27) 



This is the counterpart of Onsager's linear thermody- 
namics. The variational principle can be written in the 
form 

5A + j ^5 dydz + j dydz = 0, 

(28) 

where x D are Lagrange multipliers associated with 
the constraints given by Eqs. (P7|) . This leads to the 
relaxation equations 

^=X = -x[V + F(<7)] = -x— , (29) 



dt 



D 



+ G'{a)+^F\a) 

2y 



It is straightforward to establish that 

dydz — j -jj dydz. 



oa 



(31) 



Therefore, the relaxation equations (|29l30p satisfy A < 
if D and x a-re both positive and A> Q if D and x a-^e 
both negative. On the other hand, A = Q\EX = Y = G 
so that (^, a) is a steady state. By Lyapunov's direct 
method, we conclude that these equations can only con- 
verge towards a maximum of A (if x a-re negative) or a 
minimum of A (if D, x a-re positive). Saddle points of A 
are linearly unstable. Therefore, the relaxation equations 
(|29l30p can be used as a numerical algorithm to solve the 
optimization problem given by Eq. (1181) . 



E. Relaxation equations towards dynamical 
equilibrium 

We can introduce a set of relaxation equations that 
solve the optimization problem given by Eg. (1181) by 
adapting the general methods described in [2^ [23| ■ We 
write the relaxation equations as 



dt 



da_ 
'dt 



Y. 



(25) 



IV. BELTRAMI FLOWS 

Let us consider the minimization of energy at fixed 
helicity and angular momentum [Slf 



min{E[S„a]\H, I}. 



(32) 



The critical points of this variational principle satisfy 

5E + ii5H + aSI = 0, (33) 
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where fi and a are Lagrange multipliers. Taking the vari- 
ations of ^ and (T, we obtain 



■0 + ^(7 = 0, 

^+^i( + a = 0. 

These equations can be rearranged in the form 

1 

a = --tp, 



a 



Ji 

2/i2y /.J 



(34) 
(35) 

(36) 
(37) 



They define a steady state of the axisymmetric Euler 
equations of the form (|llll2p with f{x) = —x/jj, linear 
and g{x) ~ —ol/^ constant. In that case, the vorticity 
and the (relative) velocity are aligned everywhere 



ijJ = (u + aez X r), 



and the stationary flow is a Beltrami state. This critical 
point is a (local) minimum of energy at fixed helicity and 
angular momentum iff 

— J SS,5ijjdydz + J — dydz + ii j S£,6a dydz > 0, 

(39) 

for all perturbations 6a and (5^ that conserve helicity and 
angular momentum at first order (see Appendix [Bj. 

The variational problem given by Eq. (|32|) can be given 
several justifications: 

(i) It can be introduced in a phenomenological man- 
ner from a selective decay principle 0. Due to a small 
viscosity, or other dissipative or relaxation mechanisms, 
the energy (fragile invariant) is dissipated while helicity 
and angular momentum (robust invariants) are approx- 
imately conserved. This selective decay principle has a 
long history in physics. It first appeared in the MHD 
literature with Taylor's explanation of some behavior of 
the Zeta reversed-field pinch due to a conjectured rapid 
decay of magnetic energy relative to magnetic helicity 
[2^. This principle leads to a force- free state, i.e. a 
state whose magnetic field is proportional to its own 
curl [2^ [131. This is an analog of the Beltrami states 
that hydrodynamicists subsequently discovered in con- 
nection with axisymmetric turbulence. They are also 
related to minimum enstrophy states in 2D turbulence 
introduced by Bretherton & Haidvogel [3l| and later by 
Leith [3^ , leading to linear relationship between vorticity 
and stream function. Using the Chandrasekhar-Kendall 
eigenfunctions of the curl, these Beltrami states are easy 
to construct in both MHD and hydrodynamics. For ex- 
ample, they were used as a Galerkin basis for an exten- 
sive set of turbulent MHD computations by Shan et al. 

mp. 

(ii) In Sec. IVIl we shall propose a justification of the 
minimization problem given by Eq. p2p based on sta- 
tistical mechanics arguments. To our knowledge, this 



statistical mechanics justification has not been given be- 
fore. 

(iii) According to Eq. (|23p . the minimization problem 
given by Eq. (j32p -if it has a solution- determines a steady 
state of the axisymmetric Euler equations that is formally 
nonlinearly stable. 

Remark: the minimization problem given by Eq. (|32p 
may not have a solution, i.e. a minimum of energy 
at fixed helicity and angular momentum may not exist. 
This is the conclusion that we shall reach in [2l|. The 
absence of equilibrium state is usually associated with 
a "collapse" like the gravothermal catastrophe or the 
isothermal collapse in self-gravitating systems [35l.l36j. In 
the present context, the "collapse" is associated with the 
break-up of large scale structures and the cascade of en- 
ergy at smaller and smaller scales. The relaxation equa- 
tions associated with the minimization problem given by 
Eq. ([32]) . derived in Sec. ID 3[ may give a qualitative 



(38) idea of how the system evolves by dissipating energy 21 



However, since these equations arc purely phenomenolog- 
ical, we stress that they may not necessarily provide an 
accurate description of the true evolution of the system. 



V. STATISTICAL MECHANICS OF THE 
AXISYMMETRIC EULER-BELTRAMI SYSTEM 

A. Basic set-up 

In the previous sections, we considered steady states of 
the axisymmetric Euler equations. In a realistic situation 
where the system is forced and dissipated at small scales, 
these steady states describe the mean flow resulting from 
the balance between forcing and dissipation. However, 
there also exists fluctuations around the mean flow so 
that the velocity field is u = u -I- u', where u is the av- 
eraged velocity field and u' the fluctuations. We shall 
assume that u is axisymmetric and that the total sys- 
tem evolves while conserving the energy E = J c?r, 
the helicity H = J uTuJdr and the angular momen- 
tum I = J rugdr, but no other constraint. We assume 
that these conservation laws are selected by the proper- 
ties of forcing and dissipation, and consequently by the 
Reynolds numbers. We shall call such flows an Euler- 
Beltrami system. It was found experimentally 0, 
that the system approaches a Beltrami state when the 
Reynolds number is sufficiently large, giving support to 
the basic assumption of our theory. 

In the sequel, it will prove useful to operate 
a poloidal/toroidal decomposition such that Up = 
{ur, 0, Uz) and Ut = (0, ue, 0). In term of these fields, the 
kinetic energy density = -|- while the helicity 
density u • a; = Up • + Ut ■ uJt. For axisymmetric fields, 
j VLp-uipdv = J Uf u^t dr. To determine the distribution 
of angular momentum a and vorticity ^, we shall use a 
Mean Field Theory (MET). This method is traditionally 
very efficient in systems of high dimensionality, or with 
long-range interactions, a condition met in fluid mechan- 
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ics. In our system, we have at our disposal two privileged 
directions: the toroidal direction and the poloidal direc- 
tion. We therefore derive two different MFT procedures, 
freezing the fluctuations in one of the two directions to 
capture the fluctuations in the other direction. In each 
case, we introduce a suitable entropy, and maximize it 
under the energy, hclicity and angular momentum con- 
straints so as to obtain the Gibbs states. From these 
Gibbs states, we derive relations for the mean flow and 
for the fluctuations. The first approach, which is closely 
related to the approach in 2D turbulence will give us 
the mean field (Beltrami) and the distribution of angu- 
lar momentum (Gaussian). It will allow us to justify a 
principle of minimum energy at fixed helicity and angular 
momentum. The second approach will give us the same 
mean field and the distribution of vorticity (Gaussian). 
Fluctuations are however found to diverge in the limit of 
number of modes going to infinity, a pathology that can 
be traced back to vorticity stretching. 



B. Mean field approximation on the poloidal field: 
the distribution of angular momentum a 



1. Computations 



Let us first assume that the fluctuations are mainly in 
the toroidal direction, so that the poloidal fluctuations 
can be ignored |Up| ^ |up|. In that case, the poloidal 

held is only determined by ^ = ^ and ip, so that u' is 
made only by fluctuations of a. Let us introduce the 
density probability p(r, rj) to measure a — rj aX position 
r = {y,z). Then, the local moments of the angular mo- 
mentum are cr" = / p?7" drj. To proceed further, we need 
to introduce an entropy. Since the angular momentum 
density cr is conserved by the flow but undergoes a com- 
plicated mixing process (like the vorticity in 2D), it is 
natural to introduce the mixing entropy 



S[p] 



plnpdydzdr], 



(40) 



similar to the one introduced by Miller-Robert-Sommeria 
in 2D turbulence. We expect the entropy to increase 
during the dynamics (while the helicity, the angular mo- 
mentum and the microscopic energy are conserved) until 
the flow achieves a steady state. The functional given by 
Eq. (|40|) can also be interpreted as the neg-information 
(the opposite of the information). Maximizing this neg- 
information under given constraints is the simplest proce- 
dure we can adopt to compute the fluctuations, according 
to the information theory and its application to statisti- 
cal mechanics developed by Jaynes [37 1. 



In our approach, the conserved quantities are 
1 r 



Ef-9- ^ 



^ ^ dydz 



Ay 



dydz 



^ j ^Tpdydz + J p^dydzdri, 



H = £^adydz= / S^prjdydzdrj 



I = 



cr dydz = I prj dydzdrj. 



(41) 
(42) 
(43) 



The flrst constraint given by Eq. (|4ip will be called the 
microscopic (or flne-grained) energy because it takes into 
account the fluctuations of cr. It is different from the 
macroscopic (or coarse-grained) energy 



^ ■)/' dydz 



Ay 



dydz, 



(44) 



which ignores these fluctuations. We have E^-^- ~ 
]7jc-g- _j_ Efiuct- In our terminology, the energy will be 
called a fragile constraint because it cannot be expressed 
in terms of the coarse-grained field since cr^ 7^ ct^. While 
the microscopic energy E-^ a- jg conserved, the macro- 
scopic energy E'^ a- jg j^ot conserved and can decay. There 
is the same distinction between the fine-grained enstro- 
phy F: 
^2 — 



f-g- 



J ui'^ dr and the coarse-grained enstrophy 
/ dr in 2D turbulence [22I, 111] . On the other 
hand, the helicity given by Eq. (|^^ and the angular 
momentum given by Eq. (j43p will be called robust con- 
straints because they can be expressed in terms of the 
coarse-grained fields. We shall come back to this impor- 
tant distinction in Sec. I VII 

The most probable distribution at metaequilibrium is 
obtained by maximizing the mixing entropy S[p] at fixed 
Ef-s-, H, I and local normalization J pdrj = 1. Introduc- 
ing Lagrange multipliers, the variational principle can be 
written as 



SS - P^SE^-s- - fi^SH - a^SI 
C(r)<5 ( [ pdTj) dr = 0, 



(45) 



where is the inverse temperature and fi^ the helical 
potential (we have written these quantities with a sub- 
script <^ to recall that the fluctuations of ^ are ignored in 
the present approach). The variations on ^ imply 

/35V + MC^ = 0, (46) 
while the variations on p yield the Gibbs state 



(47) 



where the "partition function" is determined via the nor- 
malization condition 



Z(r) = / e^-^^^^^^+^^^^ci??. 



(48) 
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From Eq. (|47)) . the local average of the angular momen- 
tum is 



2y, - 
cr = - — (^j^ + a^). 



(49) 



Together with Eq. (pS)) . this equation determines a Bel- 
trami state. On the other hand /9(r, rj) , the distribution of 
the fluctuations of a, is Gaussian with centered variance 



-2 22/ 

The Gibbs state can be rewritten 

1/2 3 



V47r2/ 



(50) 



(51) 



Therefore, our statistical theory based on the conserva- 
tion of E^-^-, H and / predicts that the mean flow is 
a Beltrami state with Gaussian fluctuations of angular 
momentum. 

Note that Eq. ([50)1 means that the toroidal velocity 
fluctuations are uniform 



ue 



1 



(52) 



Therefore, can be interpreted as an inverse tempera- 
ture measuring the fluctuations of uq. These predictions 
enable the measurements of effective temperatures of tur- 
bulence through fluctuations of ue in a Beltrami flow. Be- 
cause variances are positive, f3^ is always positive, unlike 
in the 2D situation where the temperature can be neg- 
ative (in the present context, the inverse temperature is 
the equivalent of the Lagrange multiplier associated with 
the conservation of microscopic enstrophy in 2D turbu- 
lence, which is positive [l^]) . Note that Eq. ([5^ predicts 
uniformity of azimuthal velocity fluctuations which is an 
interesting prediction of our theory. This has been con- 
firmed experimentally in Q. 

On the other hand, the energy contained in the fluctu- 
ations is simply 



2. Comments 



The statistical equilibrium state given by Eqs. p6|) and 
(|49)) is of the form of Eqs. (fTTj) and (fT2|) where / is linear 
/(x) = Aa; (with A ~ — /S^/^j) and g is constant (with 
g = —a^l pL^). This means that the equilibrium state is 
a stationary solution of the axisymmetric Euler equation 
and takes the shape of a Beltrami state (see Eq. ([55]) '). 

We can also provide an interesting interpretation of 
our fluctuation relation Eq. (|52p . predicting uniformity 
of azimuthal velocity fluctuations. This equation shows 
that the azimuthal velocity fluctuations define an effec- 
tive statistical temperature l//3{. This equation may be 
regarded as formally analogous to a Fluctuation Dissipa- 
tion Relation (FDR) since it links fluctuations and tem- 
perature. These predictions enable the measurements of 
turbulence effective temperatures through fluctuations of 
ug in a Beltrami flow. As discussed previously, is al- 
ways positive. In contrast, can take positive or nega- 
tive values, depending on the helicity sign. 

The analogy between our predictions and FDRs can ac- 
tually be pushed forward. Indeed, another possible way 
to derive Eq. ([5^ is to introduce, as in classical statis- 
tical mechanics, the partition function Z describing the 
Beltrami equilibrium state in the mean fleld approxima- 
tion: 



1 S^logZ 

2 ^2 



1 5a 

tH 5l ' 



(56) 



where 5 stands for functional derivative. Formally, the 
mathematical object 6w/S^ can be seen as a response 
function. With this point of view, Eq. again reflects 
a formal analogy with FDRs since another classical way 
to write it down is to link the fluctuations of a fleld to 
its response to a perturbation. 



E 



fluct 



4y 



0-2 



1 



dydzd-q = / — dydz = —. 

Ay 2/3^ 



(53) 



Therefore, the statistical temperature ^ can also be 
interpreted as the energy (by unit volume) of the toroidal 
fluctuations. Moreover, for simple Beltrami flows with 
= 0, there is equipartition between the macroscopic 
energy in the poloidal and toroidal directions 



E/s- 



2 ^ 7 4y ^ 



with a simple connection with the helicity as 

h 



J? e.g. 



(54) 



(55) 



C. Mean field approximation in the toroidal 
direction: the distribution of vorticity ^ 



1. Spectral approach. 



We now assume that the fluctuations in the toroidal 
direction are frozen so that a — W, and that they do not 
depend on the azimuthal direction. Since the vorticity ^ 
is not conserved, we cannot in principle rigorously apply 
a statistical mechanics to the fluctuations of vorticity. 
We present here a phenomenological approach, based on 
neg-information rather than mixing entropy, and will test 
its relevance by com par ison with experimental data in 
companion papers (ol. llC|. 
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In our approach, the conserved quantities are 

^ r /• —2 



H = (T£_dydz. 



I = adydz. 



(58) 
(59) 



In the expression of the energy, we note that the fluctua- 
tions of angular momentum have been neglected so that 
cr^ = CT^ while the fluctuations of potential vorticity have 
been taken into account so that tp^ ^ tp In order to deal 
with the term -0^ that introduces a nonlocality, we shall 
develop the statistical theory in the spectral space by us- 
ing an approach similar to that developed by Kraichnan 
[39| and Salmon et al. [S^] in 2D turbulence. Let us 
first note that the field cj) = ip/r satisfies the differential 
equation 



1 



(60) 



To solve the problem, we decompose the fields onto the 
eigenfunctions (j)mn of C defined by 



rnn i 



B;an4>mn, (61) 



with 4>mn = on the boundary. Taking the origin of the 
z axis at the bottom of the domain, the eigenfunctions 
are given by the Hankel-Fourier modes 



where jim is the m*'' zero of Bessel function Ji. The 
mode (pmn corresponds to m cells in the radial direction 
and n cells in the vertical direction. The corresponding 
eigenvalues are 



Jim \ , Cmry 



R \2hJ 



(63) 



The eigenfunctions are orthogonal with respect to the 
scalar product 

ifg) = / / rdrdzfg, (64) 
Jo Jo 

so that {(I)mn4>m'n') = ^mm'&nn' whcre Jj-, is the Kro- 
necker symbol. We now decompose the fields on these 
eigenmodes writing 



t 
r 



EE* 

m— 1 n— 



m—1 n—0 

cr 

U0 



(65) 
(66) 
(67) 



where we have restricted the sum over finite number of 
modes, so as to respect Hamiltonian condition for the 
Euler equation. Moreover, the finite number of modes 
implies a coarse graining of the solution, that is desirable 
to reach a stationary state. With this decomposition, 
using Eq. (|60p . we have by construction 



Bmn^mn- 



(68) 



Inserting the decomposition given by Eqs. (|65II67I) in the 
expression of the energy and helicity, and using the or- 
thogonality condition given by Eq. ([M)) and the Parseval 
identities, we obtain 



E =-^( 



mn I 2 



mn 



(69) 
(70) 
(71) 



where the brackets denote a domain average. To apply 
the statistical theory, we introduce the density probabil- 
ity Pmn{i^) of measuring the value ujmn = v of the vortic- 
ity in the mode {m,n). We can rewrite the constraints 
in the form 

^ = ^ E / + '^™"(^) '^"^ (^2) 

H = ^ j w,„„u,„„p„i„(i/)di/, (73) 

mn 

I = ^U„in{r(l)„in)- (74) 
mn 

2. The multi-modes case 

We first consider the situation where the energy is 
spread all over the different wavenumbers {n,p). In that 
case, there is no reason to specialize a special mode and 
we can use the neg-information in spectral space to define 
the entropy 



5 = -^ j Pmn{j^) ln/9m„(z^) dv. 



(75) 



The Gibbs state is obtained by maximizing S at fixed 
energy, helicity, angular momentum and normalization 
/ Pmn dv = 1. We write the variational principle as 



- ^ Xmn^ i / Pmndv j = 0, 

mn / 



(76) 



m—l n—0 



where p^ and are Lagrange multipliers (the sub- 
script (T recalls that the fluctuations of angular momen- 
tum have been neglected). The variations on Umn yield 

/3crWmn + Pa^^mn + acr{r(l>mn) = 0. (77) 
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Multiplying by (f>mm summing on the modes and using 

= Y.mn{'^^ran)4'mn, we obtain 



— h llat, + = 0. 

The variations on yield the Gibbs state 



P7n7i{v) 



l-MirMmTif 



(78) 



(79) 



where Zmn is a factor ensuring the local normalization 
condition. The distribution of the fluctuations of ujmn is 
therefore Gaussian with 



Po-— H Mo- tin 



52 

mn 



0, 



d2 

-2 ran 



(80) 

(81) 



According to Eq. (jSip . the mean fluctuating energy per 
mode number (w^„ — cJ^)^/i3^„„ is constant: we have 
equipartition of energy for the poloidal fluctuations. This 
result is a classical outcome of equilibrium statistical me- 
chanics. Multiplying Eq. (j80p by 4>mn, summing on the 
modes and using Eq. (|68p . we obtain 



(82) 



Equations ([75]) and ([5^ show that the mean flow associ- 
ated with the statistical equilibrium state is a stationary 
solution of the axisymmetric Euler equations correspond- 
ing to a Beltrami state. One can also deduce from Eq. 
([8T|l that the total volumic energy contained in the fluc- 
tuations is 



E 



fluct 



E 



{v - w„„)2 Ntot 



2B2 



(83) 



where iVtot is the total number of modes. The linear 
divergence with Ntot comes from energy equipartition. 
Therefore, the statistical temperature in this case is pro- 
portional to the poloidal energy of fluctuations. This is 
therefore analog to the previous mean field case. 

It is also interesting to compute the azimuthal vorticity 
fluctuations. They are given by 



2 i2 
mn r^mn i 



(84) 



since uimn and ujm'n' are independent if (m, n) ^ (m', n'). 
For comparison with real turbulent data fields, we have 
constructed synthetic instantaneous fields obeying Eqs. 
([TS]). ([7^ and to study some of their properties. An 
example is shown in Fig. [1] The average velocity field is 
obtained as follows: we choose a given number of modes 
to represent the field (A^,„, Nn)- Then, we get the coeffi- 
cient of the Hankel-decomposition Umn by a least-square 
fit to ug of an actual mean turbulent field (here, the ve- 
locity field obtained by counter-rotation at F = &Hz of 



TM73 impellers-this field is described in Q). This field is 
shown in Fig. [T] Once the Umn are obtained, we get the 
coefficient of the Hankel-decomposition of the mean 
stream function from a least-square fit to the poloidal 
experimental velocity fields u^. and using Eqs. 
Such a field, shown in Fig. [1] obeys the relation ([5^ 
with B = —fia/pa = —3.6, as shown in Fig. [2] We 
then obtain the mean azimuthal vorticity LOg thanks to 
((60)) . It is shown in Fig. [1] We also use the modes of 
the decomposition to compute the theoretical variance 



following Eq. 



, shown in Fig. [21 Finally, we 



compute an instantaneous fluctuation fleld of azimuthal 
vorticity by drawing for each mode a realization of the 
Gaussian distribution Eq. ([7^ and reconstructing the 
field through Eq. ([55)) . The results are provided in Fig. 
[T] Wc see that the theoretical variance w^ms not inde- 
pendent of r and z but its dependence in z is weak. In 
the radial direction, it oscillates mildly around a value 
■ 



(85) 



Noteworthy, is much larger than the value B"^ / P^r 
obtained in the one mode case (see Eq. (jlOOp of Sec. 
IV C 3[) . This is reminiscent of what has been observed in 
real turbulent data fields It is therefore interesting 
to study further the dependence of Q.00 with respect to 
the number of modes of the problem. We have found 
empirically that this value behaves like 



^00 — Ntot- 



(B^ 



Pa 



where {B"^) is a mean Beltrami factor, defined as: 



1 



Nt, 



2 

mn ■ 



(86) 



(87) 



This number depends on the set of modes (values of n 
and m) considered. For example, in an isotropic situation 
when one sums only over modes such that n — m, i?^„ ^ 



and (^2) 



Nt-lN^ 



Nl^t' (here Nt, 



N'^). In an anisotropic situation such as A^,„ = I7 ^ 
n2 and (^2) ~ N^^IN^ ^ Nit (here Ntot = N^N^ = 
Nn). 

Integrating Eq. ([84)) over the volume, we get the az- 
imuthal enstrophy fluctuations 



fluct 



(a;2 _ ujg^)dydz = ^ ^ 

Per 



Nt, 



Pa 



IB' 



Since {B'^) grows algebraically with Ntot (see above), the 
enstrophy fluctuations therefore diverge as Ntot x -^tot 
when the number of modes becomes inflnite. The lin- 
ear part of the divergence is comparable with the linear 
divergence obtained for the variance of energy fluctua- 
tions, and can be thought to be an outcome of equilib- 
rium statistical mechanics. Note that in real turbulent 
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FIG. 1: (Color online) Example of synthetic instantaneous 
fields in the multi-mode case, at —fialiia ~ —3.6 and Pa ~ 2: 
(a) projection of {lu, ue, IZI) over a vertical meridional plane. 
The toroidal velocity is coded in color, the poloidal velocity 
is denoted with arrows; (b) projection of (p over a vertical 
meridional plane; (c) projection of Uo over a vertical merid- 
ional plane; (d) projection of uj'g over a vertical meridional 
plane. The field is constructed with Nm = 2 radial modes 
and Nm = 2 axial modes. The amplitude of the velocity 
modes results from a least square fit to the average velocity 
field described in with a forcing using TM73 turbines in 
contra-rotation at F = GHz over eigenmodes given by Eqs. 
(|61l) . The instantaneous value is obtained by drawing one 
realization of random vorticity fluctuations according to the 
theoretical Gibbs distribution given by Eq. (|79p . 



flows, -out-of-equilibrium solutions with non-zero energy 
flux-, the energy of fluctuations by mode Efiuct/Ntot re- 
mains finite but the enstrophy of fluctuations by mode 
^fiuct/Ntot ^ ^uit diverges algebraically with the num- 
ber of modes in 3D, while the divergence is much milder 
(logarithmic) in 2D. This difference can be seen as the 
signature of 3D vortex stretching, that is captured by 
our model. Indeed, vortex stretching induces a trans- 
fer of vorticity at continuously decreasing scales, thereby 
leading to an enstrophy divergence. 

Note that both the enstrophy fluctuations and the en- 
ergy fluctuations provide an estimate of the statistical 
temperature. The comparison of temperature in between 
the two measurements actually provides a measure of the 
number of modes in the system since energy fluctuations 
are proportional to Ntot and enstrophy fluctuations be- 
have roughly like iV^^f for an isotropic situation. This 
will be further discussed in Sec. IV C 41 



3. The one mode case 



4 

a" 
-4 

(a)-l 




-0.5 0.5 




FIG. 2: (a) (Color online) Toroidal velocity ue as a function 
of the stream function 4>. The line is a linear fit, with slope 
B — —3.6. (b) Spatial variation of the theoretical vorticity 
variance as a function of r only. This variance has been ob- 
tained using the modes of the Hankel decomposition of the 
velocity field and Eq. (|84l) . The dispersion in vertical di- 
rection corresponds to a dispersion of variance at a given r 
with varying z. The horizontal line is the empirical value 
Ntot{B'^) / Per ■ The conditions are the same as in Fig. [T] 



of entropy). It is interesting to consider the opposite 
case, where all the energy is concentrated in one mode 
(toc'^c)- The results pertaining to this case have been 
discussed in We present here the corresponding de- 
tailed computation. In such a case, it is easy to check 
that 



(89) 



where we have noted B = Bm^n^ for brevity. Fur- 
thermore, the term Ky can always be introduced since 
it is in the kernel of A. Since the energy is concen- 
trated in one mode in the spectral space, we cannot 
use the neg-information anymore in that space to de- 
termine the Gibbs distribution. However, by the Lin- 
quist theorem, such a peaked probability in the spectral 
space corresponds to a spread probability in the physical 
space. We therefore turn back to the density probabil- 
ity p(r, v) to measure ^ = at position r, introduce the 
neg-information 



S[p] 



p In p dydzdv, 



(90) 



in the physical space, and consider the maximization of 
S[p] at fixed E, H and normalization J pdv — 1. Using 
Eq. ([Sn]) . the constraints can be written 



— dydz 
Ay 



V 



= I '^P^ dydzdv + — j ypvdydzdv+ j — dydz, 



Ay 



(91) 



The previous Gibbs distributions do not couple modes 
with different wavenumbers. They are the analogs of 
the Boltzmann laws found in simple quasi-geostrophic 
models [40| in the case where energy is spread evenly 
over all modes (an implicit assumption behind our choice 



H 



a ^ dydz = / apiy dydzdv, 



I ^ a dydz. 



(92) 



(93) 
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We write the variational principle as 

SS — /ScrSE — fiaSH — a^SI 
x(r)(5 I / pdiy] dr = 0. 



The variations on a imply 

-X h MctC + a<T = 0, 

and the variations on p yield the Gibbs state 

z 



(94) 



(95) 



(96) 



where Z is the normalization factor. The distribution is 
therefore Gaussian. Its first two moments are 



B 



^2 ^ e-^ = 



2P-a O- B^K 



The Gibbs state can be rewritten 



a \ 1/2 



7rB2 



(97) 
(98) 

(99) 



One sees that relations (|97l) and (|95|) can be satisfied si- 
multaneously only if B'^ = {fi^j [laY and — ^oiap,al P^. 
This fixes the wavcnumbcrs {rric, ric) of the mode in which 
energy is accumulated and provides a physical interpreta- 
tion of the ratio Pa/fJ-a- The relations (|97|) and ([89|) show 
that the mean flow associated with the statistical equi- 
librium state is a stationary solution of the axisymmetric 
Euler equation corresponding to a Beltrami state. More- 
over, Eq. ((98|) shows that the fluctuations of vorticity are 
uniform and scale as 



9 2 

- UJg 



B^ 



(100) 



Integrating over the volume, we find that the azimuthal 
enstropliy fluctuations are 



Per 



(101) 



We may also compute the energy contained in the fluc- 
tuations as 



E 



fluct 



S2 



pdydzdrj = 



l^^,dydz^^^. 



(102) 



Therefore, the statistical temperature can be sim- 
ply interpreted here as the volumic energy of the 
poloidal fluctuations. Note that in such a simple case, 
^ fluct I E fluct = 1 consistent with the multi-mode re- 
sult since B^ = (B'^) in the "one mode case". Moreover, 



for simple Beltrami flows with acr = 0, there is equiparti- 
tion between the energy of the mean flow in the poloidal 
Ep and toroidal direction Et: 



1 /"-r^ 



with a simple connection with the mean helicity as 



(104) 



This case is therefore the analog to the mean field toroidal 
case, with uniformity of fluctuations and simple connec- 
tion with helicity. 



i. A note on vorticity fluctuations 

The link of the present results with experiments has 
been partially discussed in 0, Q , with focus on the mean 
field toroidal case and the one-mode mean field poloidal 
case. A detailed comparison is provided in a compan- 
ion paper [l3|. It is however interesting to come back to 
one puzzling result of @ to give it a new interpretation 
in the present context. It has indeed been found that 
in a turbulent counter-rotating von Karman flow both 
the velocity and vorticity fluctuations are approximately 
uniform over the box. When interpreted in the context 
of the mean field toroidal case, and the one-mode mean 
field poloidal case, the value of this constant provides 
an estimate of the two inverse temperatures /3j and I3„ 
through Eqs. ((52|) and pOOp . Experimentally, one finds 
/3j ^ Pa, with a ratio (3^ /Pa ranging from 8 to 17 in 
different forcing configurations. In the restricted context 
discussed in Q, this difference is puzzling, and points 
towards the existence of two different temperatures. If 
however one considers a wider context, in which the vor- 
ticity fluctuations are considered to span several modes, 
the experimental measurements allow for another inter- 
pretation. Indeed, using our least-square fitting to exper- 
imental data with varying Nn and Nm with iV„ ~ N^a 
(isotropic case), we found that — uJe'^ ~ ^tot^f~^ in- 
stead of the value 3.6^/ p^ predicted in the one-mode case 
Ntot ~ 1 (recall that B = —3.6 in the selected experimen- 
tal data). Assuming P^ = Pa- instead of Ntot = 1, we can 
then use the experimental measurements to infer Ntot- 
One finds a value ranging from Ntot = 4 to Ntot = 6. 
Such a small number is intriguing because the usual be- 
lief is that turbulent flows are characterized by a very 
large number of degrees of freedom. There are however 
other indications {e.g. in the dynamo context [4l| ) that 
turbulent flows can be described using tools adapted from 
dynamical systems, as if the effective number of degrees 
of freedom were indeed small. 
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VI. CONNECTION BETWEEN DIFFERENT 
VARIATIONAL PRINCIPLES 

In this section, we make the connection between dif- 
ferent variational principles that characterize the equilib- 
rium states. For each principle, we propose a relaxation 
equation (see Appendix |D| that can be used as a nu- 
merical algorithm to solve the corresponding variational 
problem. These relaxation equations can also provide 
an effective description of the relaxation of the system 
towards the equilibrium state. They will be solved nu- 
merically in [1^. We finally justify through statistical 
mechanics the phenomenological principle according to 
which: "the mean flow should minimize the macroscopic 
energy at fixed helicity and angular momentum" . In this 
section, we ignore the fluctuations of vorticity and exclu- 
sively consider the mean field theory of the poloidal field 
developed in Sec. IV Bl 



with 



2y - 



a2 ~ — 



(113) 
(114) 



As already indicated, the last relation shows that the 
temperature must be positive. Furthermore, a critical 
point of (fTOS]) is an entropy maximum at fixed E-^'^' , 77, 
/ and normalization iff 



dydzdrj 



S^StP dydz - fi / S^Sadydz < 0, (115) 



for all perturbations Sp and 6^ that conserve energy, he- 
licity, angular momentum and normalization at first or- 
der . 



A. The basic variational principle 

The basic maximization problem that we have to solve 



is: 



with 



and 



max{5'[p] I E^-a-, H, /, / pdr] = 1}, (105) 



S[p] 



pin p dydzdrj, 



(106) 



Ef-s- ^ - J ^tpdydz + J — dydz 

-,2 



= ^ J^ipdydz + J p^dydzdr], (107) 
H = J S^lf dydz = j S^prj dydzdrj, (108) 
I = Wdydz^ / prjdydzdrj. (109) 



The critical points arc determined by the variational 
principle 

5S - fi5E^-3- - fi5H - a5I 
C{r)S ( I pdij) dydz = 0. (110) 



(111) 



The variations on ^ imply 

while the variations on p yield the Gibbs state 

/ o \ 1/2 

p{r,rj) = i^-£-j e-i(''--)^ 



(112) 



B. An equivalent but simpler variational principle 



The maximization problem (|105p is difficult to solve 
because the stability condition ()115|) is expressed in terms 
of the distribution p{r,r]). We shall here introduce an 
equivalent but simpler maximization problem by "pro- 
jecting" the distribution on a smaller subspace. To solve 
the maximization problem (|105p . we can proceed in two 
steps [H- 

(i) First step: we first maximize S at fixed E-^'^', H, I, 
J pdr] = 1 and a{r) ~ J prjdrj and ^(r). Since the spec- 
ification of CT(r) and ^(r) determines / ip£,dr, H and /, 

2 

this is equivalent to maximizing S at fixed } dydzdrj, 

J pdrj = 1 and a{r) — J prjdrj. Writing the variational 
problem as 



dS — pS [ / p— dydzdrj 

4y 



\{r)5 I prjdrj ) dydz 



C{r)6 pdrj] dydz = 0, (116) 



we obtain 



(117) 



and we check that it is a global entropy maximum with 
the previous constraints since S^S ^ — J dydzdrj < 
(the constraints arc linear in p so their second variations 
vanish) . We also note that the centered local variance of 
this distribution pi is 



(72 = o"-^ — cr^ 



2y 

/3' 



(118) 



implying /3 > 0. Using the optimal distribution given by 
Eq. (|117|) . we can now express the functional in terms 
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of a" and /3 writing S = S[pi] and = Ef-9-[pi] 

After straightforward calculations, we obtain 



5=--ln/3, 



(119) 



Ef-s- ^IJ ^^dydz + ^ + dydz, (120) 



H = £,(7 dydz, 



a dydz, 



(121) 



(122) 



up to some constant terms. Note that /3 is determined 
by the energy constraint given by Eq. (fT20| leading to 



1 



(123) 



where iJ'^f ■ is the macroscopic energy defined by Eq. 

This relation can be used to express the entropy 
Eq. pi9p in terms of ^ and a alone. 

(ii) Second step: we now have to solve the maximiza- 
tion problem 



ma_x{5[a,e]|i?^-^-,i?, /}, 



(124) 



with 



5 = i In {2Ef-3- - j dydz - dydz ) , (125) 



H 



i a dydz, 



a dydz. 



(126) 



(127) 



(iii) Conclusion: Finally, the solution of (jlOSp is given 
by Eq. (|117p where a is solution of (jl24p . Therefore, 
(jlOSp and (|124p are equivalent but (|124|) is easier to solve 
because it is expressed in terms of a and ^ while (|105|) is 
expressed in terms of p and ^. 

Up to second order, the variations of entropy given by 
Eq. ([T^ are 



AS* 



/3 2 / iljSidydz+ / 5tp5£,dydz 



mi 

22/ 



dydz 



u6c 



V 



dydz 



abc 



-/32( / ^didydz+ I ^dydz] , (128) 



where /3 is given by Eq. (|123|) . The critical points of 
(|124p are determined by the variational problem 



5S - fiSH - aSI = 0. 



The variations on ^ yield 



l3il; + pW = 0, 



and the variations on a yield 



(129) 



(130) 



(131) 



This returns the equations (|llip . (|113p for the mean flow. 
Together with Eq. (|117p . we recover the Gibbs state 
given by Eq. ()112p . Considering now the second varia- 
tions of entropy given by Eq. (|128p , we find that a critical 
point of (jl24p is a maximum of S at fixed microscopic en- 
ergy, helicity and angular momentum iff 



1 



/3 / S^jjSi dydz 



2y 



dydz ] ~ p I S^Sa dydz 



^pSi dydz 



dydz < 0, (132) 



for all perturbations and Sa that conserve helicity and 
angular momentum at first order (the conservation of mi- 
croscopic energy has been automatically taken into ac- 
count in our formulation). The stability criterion (jl32p 
is equivalent to Eq. pisp but it is much simpler because 
it depends only on the perturbations and SW instead 
of the perturbations 6p of the full distribution of angular 
momentum. In fact, the stability condition p32p can be 
further simplified. Indeed, using Eqs. (|130p and (|13ip . 
we find that the last term in parenthesis can be written 

ipSi dydz + / dydz = 

J 2y 

J {w6i + lSa) dydz-^ J 5a dydz, (133) 

and it vanishes since the helicity and the angular mo- 
mentum are conserved at first order so that SH = 
J {a6i + i6a) dydz ~ and 61 = J 6 a dydz ~ 0. There- 
fore, a critical point of (|124p is a maximum of entropy at 
fixed microscopic energy, helicity and angular momentum 
iff 



/ Si^Sidydz 



22/ 



■ dydi 



p I S^Sa dydz < 0, 



(134) 



for all perturbations 6^ and 6a that conserve helicity and 
angular momentum at first order. In fact, this stability 
condition can be obtained more rapidly if we remark that 
the maximization problem (|124p is equivalent to the min- 
imization of the macroscopic energy at fixed helicity and 
angular momentum (see Sec. IVI Cp . 
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C. Equivalence with the minimum energy principle 

Since ln(a;) is a monotonically increasing function, it is 
clear that the maximization problem (11241) is equivalent 
to 



mm{E''-3 [a,^]\H, I}, 



(135) 



with 



^c.g. f ^^dydz+ f — dydz, (136) 
2 J J 4y 



H 



^ a dydz, 



I = J a dydz. 
We have the equivalence 



(137) 



(138) 



(139) 



Therefore, the maximization of entropy at fixed micro- 
scopic energy, helicity and angular momentum is equiv- 
alent to the minimization of macroscopic energy at fixed 
helicity and angular momentum. The solution of (|105l) is 
given by Eq. ()117p where a is solution of p35p . There- 
fore (|105p and (jl35p are equivalent but (|135p is easier 
to solve because it is expressed in terms of W instead of 
p. Our approach therefore provides a justification of the 
minimum energy principle in terms of statistical mechan- 
ics. Note that, according to ([25]) . the principle (|135p also 
assures that the mean flow associated with the statistical 
equilibrium state is nonlinearly dynamically stable with 
respect to the axisymmetric Euler equations. 

The critical points of (|135p are given by the variational 
problem 

SE"-^- + fiSH + aSI = 0. (140) 
The variations on ^ yield 

ij + fia^O, (141) 
and the variations on a yield 

a^~2y{fil+a). (142) 



This returns Eqs. ()111|) and (jll3|l for the mean flow (up 
to a trivial redefinition of fj, and a). Together with Eq. 
(|117p . we recover the Gibbs state given by Eq. (|112p . On 
the other hand, this state is a minimum of at fixed 
H and / iff 



i / 6'il;5£,dydz 



{6af 



Ay 

fi I SS^SWdydz > 



dydz 



(143) 



for all perturbations and Sa that conserve helicity and 
angular momentum at first order. This is equivalent to 
the criterion given by Eq. (|134p as it should. 

We have thus shown the equivalence between the max- 
imization of Boltzmann entropy at fixed helicity, angular 
momentum and fine-grained energy with the minimiza- 
tion of coarse-grained energy at fixed helicity and angu- 
lar momentum. This equivalence has been shown here 
for global maximization. In Appendix [Cl we prove the 
equivalence for local maximization by showing that the 
stability criteria (|115p and (|143p are equivalent. 



D. Equivalence with the canonical ensemble 



The basic maximization problem ()105p is associated 
with the microcanonical ensemble since the energy 
is fixed. We could also introduce a canonical ensemble 
where the inverse temperature 13 is fixed by making a 
Legendre transform J = S — (3Ef-s- of the entropy with 
respect to the energy (s^ . The corresponding maximiza- 
tion problem is 



max{ J[p] \H, I, / pdrj = 1}. 



(144) 



A solution of (jl44p is always a solution of the more con- 
strained dual problem (jl05p but the reciprocal is wrong 
in case of ensembles inequivalence. In the present case, 
however, we shall show that the microcanonical ensem- 
ble (|105p and the canonical ensemble (jl44p are equiva- 
lent. This is because the fluctuations of the energy are 
quadratic. 

To solve the maximization problem (|144p we can pro- 
ceed in two steps. We first maximize J at fixed I, 
J pdrj = 1 and a{r) = / p7]dri and ^(r). This is equiv- 
alent to maximizing J = S ~ (3 J p^ dydzdrj at fixed 
J pdrj = 1 and ^(r) = J prj drj. This leads to the optimal 
distribution (|117p where /? is now fixed. This is clearly 
the global maximum of J with the previous constraints. 
Using this optimal distribution, we can now express the 
free energy in terms of ^ and W by writing J[^, a] ~ ^[pi]- 
After straightforward calculations, we obtain 



J = -PE"" 



(145) 



up to some constant terms (recall that /? is a fixed param- 
eter in the present "canonical" situation). In the second 
step, we have to solve the maximization problem 



max{J[CT,e] /}. 



(146) 



Finally, the solution of (|144p is given by (|117p where 
is determined by (|146p . Therefore, the canonical varia- 
tional principle (|144p is equivalent to (|146p . On the other 
hand, since /3 > 0, the maximization problem (jl46p with 
Eq. (|145p is equivalent to (|135p . Since we have proven 
previously that (|135p is equivalent to the microcanonical 
variational principle (|105p . we conclude that the micro- 
canonical and canonical ensembles are equivalent. 
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VII. CONCLUSION 

In the present paper, we have constructed a simph- 
fied thermodynamic approach of the axisymmetric Eu- 
ler equations so as to describe its statistical equihbrium 
states. This predicts the mean field at metaequilibrium 
and the fluctuations around it. We have considered two 
mean field theories. In the first one, we have ignored the 
fluctuations of vorticity. In that case, we have found that 
the fluctuations of angular momentum are Gaussian and 
that the mean flow is in a Beltrami state. Furthermore, 
we have proven that the maximization of entropy at fixed 
helicity, angular momentum and microscopic energy is 
equivalent to the minimization of macroscopic energy at 
fixed helicity and angular momentum. This provides a 
justification of this selective decay principle from statis- 
tical mechanics. These results are very similar to the 
case of 2D turbulence if we make the analogy between 
the angular momentum (axisymmetric) and the vorticity 
(2D). Indeed, in the simplified statistical approach of the 
2D Euler equations developed in [22|, the fluctuations of 
vorticity are Gaussian and the mean flow is characterized 
by a linear Zu — tp relationship. Furthermore, it has been 
proven that the maximization of entropy at fixed energy, 
circulation and microscopic enstrophy is equivalent to the 
minimization of macroscopic enstrophy at fixed energy 
and circulation. This provides a justification of the min- 
imum enstrophy principle from statistical mechanics. In 
the second mean field theory, we have ignored the fluc- 
tuations of angular momentum. In that case, we have 
found again that the fluctuations of potential vorticity 
are Gaussian and that the mean flow is in a Beltrami 
state. We have also observed an interesting signature of 
vorticity stretching, via divergency of the variance of the 
vorticity fluctuations with increasing number of degrees 
of freedom. Overall, the variance of fluctuations provides 
a measure of the number of degrees of freedom and of 
the statistical temperature(s) of turbulence, and allows 
to check Fluctuation-Dissipation Relations (FDR). 

The question is whether these results are applicable to 
a laboratory flow such as von Karman flow. On the one 
hand, several basic hypotheses are not satisfied in the VK 
flow: it is a dissipative, forced flow, and instantaneous ve- 
locity fields are not axisymmetric. On the other hand, 
one observes that in a stationary state, dissipation and 
forcing balance globally, and may be neglected locally, 
at least within an inertial range of scales, and the mean 
flow is axisymmetric. This motivated experimental tests 
of the equilibrium and fluctuations relations, reported in 
in a large Reynolds number von Karman flow. It has 
been found that the observed stationary states are well 
described by the equilibrium states of the Euler equa- 
tions, and that both the velocity and vorticity fluctu- 
ations are approximately uniform over the box. These 
fluctuations depend on three unknowns: the mean fleld 
statistical temperatures l//?^ and l//3cr, and the effective 
number of degrees of freedom Ntot- One can therefore 
use the experimental measurements to estimate these pa- 



rameters under a supplementary assumption. Assuming 
Ntot = 1, one flnds that Pa and /3{ differ by an order of 
magnitude 0], an intriguing result that may come from 
the out-of-equilibrium character of the turbulence. On 
the other hand, assuming /S^ = P^ one flnds that Ntot is 
of the order of 4 to 6, a rather small number for a turbu- 
lent flow. This may be due to strong correlations within 
the flow, that effectively reduces the number of degrees 
of freedom. 

Altogether, these results are an indication that ther- 
modynamics of Euler axisymmetric flows can bring new 
interesting information about real flows. Indeed, our sta- 
tistical theory of axisymmetric flows that can account for 
certain experimental results reported in 0, • In a forth- 
coming communication [2l| we shall explain from this 
approach a turbulent bifurcation that has been observed 
in a von Karman flow psj . Despite its good agreement 
with experiments, several criticisms can be made to our 
approach. For example, in contrast with what is hypothe- 
sized or derived in the paper, it is indeed usually believed 
that the statistics of fluctuations in turbulence are essen- 
tially non-Gaussian, that the mean-fleld approximation 
is not satisfactory, and that the forcing and dissipation 
are important in the process. We remark that these prop- 
erties have been observed for homogeneous and isotropic 
(three dimensional) turbulence. In our 2.5D situation, 
these general results may not be correct anymore, since 
the dominant dynamical processes in both systems are 
probably different, in particular because of the quasi-2D 
nature of axisymmetric turbulence. In our case, it is not 
unlikely that the distribution of angular momentum is 
Gaussian (or close to Gaussian), like the distribution of 
the velocity components in 3D turbulence. However, the 
distribution of azimuthal vorticity, which is a derivative, 
may not be Gaussian because of intermittency. Unfortu- 
nately, we are not able to check these predictions exper- 
imentally due to a lack of statistics. 

We have applied the maximum entropy principle to 
the Euler equation. Alternatively, Adzhemyan & Nal- 
imov [4^, [131 have applied this principle to the stochastic 
Navier-Stokes equation. The renormalization group ap- 
proach was used instead of the mean-field approximation, 
a non-Gaussian distribution was obtained and the Kol- 
mogorov spectrum was derived for the inertial range. It 
would be interesting to extend their approach to our 2.5 
situation to see what it predicts and compare with our 
results. 

This work was supported by European Contract 
WALLTURB. 



Appendix A: Generalization and link with other 
results 

In the present paper, we have considered a restricted 
class of flows for which the only invariants are E, H = Hi 
and / = /i . The complete generalization to arbitrary in- 
variants of axisymmetric Euler equations (E, Hp, Iq) 
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remains an unsolved problem. There are however some 
special cases where we can perform the maximization 
problem and find equilibrium distributions. We list these 
cases below, and make connection with previous results. 



Since 172 > 0, the generalized entropy is proportional to 
minus l!^'^'- The corresponding critical points 

= (AlO) 



1. Conservation of E"'^', H = Hi,V = Ho, I = h and 

rf-9' 

'-n>l 

Leprovost et al. Q have considered the maximization 
problem 

max{S[p] I H, r, /, f pdr, = 1}, (Al) 



where S[p] is the mixing entropy and I^y'^ 

J prf^ dydzdrj. If we make a Legendre transform of the 
entropy with respect to the fragile constraints, we obtain 
the reduced maximization problem 



with 



max{S^[p] I E^-a-, H, T, /, / pdi] = 1}, 
SM-S[p]-J2arJl-^-. 



n>l 



Explicitly 



pin 



xiv) 



dydzdrj, 



(A2) 



(A3) 



(A4) 



where x('7) = e 5i;„>i "n'; _ Proceeding as in [42j, we can 
show that the maximization problem \A2\ is equivalent 
to 



mSiyi{S[a]\E''-9-, H, T, /}, 
where S[a] is the generalized entropy 



(A5) 



S = -J C{a)dr, C{a) = - j\{\nx)']-\-x)dx. 

(A6) 

The critical points of (|A5|) are given by 



-pa - 7, 



2y 



(A7) 



(A8) 



The solutions of (|A2[) and (|A5[) are always solutions of 
(|Aip but the reciprocal is wrong in case of ensemble in- 
equivalence. Thus (|X2l) If we con- 
sider the particular case where = for n 7^ 2 and 
a2 ^ 0, then, proceeding as in [i^l, we find that the fluc- 
tuations of angular momentum are Gaussian with vari- 
ance (72 = l/(2a2) and that the generalized entropy is 



2(72 



dv = —I^-'-. 

2(72 



(A9) 



0-2 2?/ 



(All) 



are steady states of the axisymmetric Euler equations 
corresponding to /('0) linear and g{^) linear but not con- 
stant. Therefore, the statistical approach of Leprovost et 
al. based on (|A1|) does not lead to Beltrami states 
(corresponding to / linear and g constant) contrary to 
the statistical approach developed in the present paper. 
Finally, if we consider the maximization problem 

ma_x{5[p] I E- a-, H, T, /, / pdi^ - 1}, (A12) 

where only the quadratic integral j/'^' is conserved 
among the set of fragile constraints, and proceed as in 
[2^ . we find that (|A1[) is equivalent to 



mm{I^-^ [a]\E^-s-,H, T, /}. 



(A13) 



Since (jAlSp is equivalent to (|A5|) with (jA9[) . hence to 
()A2p . we also conclude that, in the specific case where 



only Jj'^' is conserved among the set of fragile con- 
straints, (jA12[) is equivalent to (jA2[) with = for 
n 7^ 2. 



2. Conservation of H = Hi, T = Ho and / = /i 

We consider the maximization problem 

max{5[p] I -B^-^'-, iJ, T, /, / pdry = 1}, (A14) 



The variations over ^ give 

= -/iff - 7, 



(A15) 



and the variations over p give the exponential distribu- 
tion 



lp-(A*^j+K+«)') 
Z 



(A16) 



Since this distribution is not normalizable, we must im- 
pose some bounds on the angular momentum and we 
shall assume — A < (7 < A (symmetric). In that case, we 
have 



7f= XL 



-A(/3— + a 



where 



L{x) = tan ^{x) 



(A17) 



(A18) 
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is the Langcvin function [4^. We see that is a 

steady state of the axisymmetric Euler equations. If we 
consider the maximization problem 

max{S[p] I E'^-s-, H, T, I, [ pdv = 1}, (A19) 

we find symmetric results but, in that case, (^,0") is not 
a steady state of the axisymmetric Euler equations. 



3. Conservation of fi^ ^ , H = Hi and I = h 

In the present paper, we have considered the maxi- 
mization problem 



max{S[p] I H, I, [ pdi] = 1}. 



(A20) 



The variations over ^ give 

Pij + pW=0 (A21) 
and the variations over p give the Gaussian distribution 



1 

p = — e If 



(A22) 



We have 



^=-J{^4 + a), a,^^. (A23) 

In that case is a steady state of the Euler equa- 

tions with / linear and g constant (Beltrami state). The 
stream function can be expressed in terms of Bcssel func- 
tions. Furthermore, we have shown that (jA20p is equiv- 
alent to 



mm{E''<^ [W]\H, I}. 



(A24) 



4. General case: conservation of E^'^' , H = Hi, I, 

ljf-9- rf-a- 



Let us consider the problem 

max{5[p] I Ef-s-, /, H^^^ 1^% [ pdrj = 1}, 

(A25) 

generalizing the one studied in the present paper. The 
variations over ^ give 

P^p + pa + Y^ Pn^ = 0, (A26) 

)i>i 

and the variations over p give 

p = i_g-E„>i""')"g-E„>iMn5';"g-T^g-(Ml'+a)')^ 

^ (A27) 
However, it is difficult to be more explicit. Therefore, we 
shall consider simpler problems. 



5. Conservation of E^'^', H = Hi, I = h, and /| 
We consider the maximization problem 

max{S'[/9] I Ef-3-, H, I, I^^- , [ pdrj = 1}. (A28) 
p,? J 

The variations over ^ give 

Pi' + pa = 0, (A29) 

and the variations over p give the Gaussian distribution 



1 ('1-^)^ 

g 2<T2 ^ 



with 



pS, + a 

"2(a2 + |)' ^"'"2(^2 + ^, 



(A30) 



(A31) 



In that case (^,5") is a steady state of the Euler equa- 
tions with / and g linear. The stream function can be 
expressed in terms of Whittaker functions Q . If we make 
a Legendre transform of the entropy with respect to the 
fragile constraints, we obtain the reduced maximization 
problem 

max{5[p] \H,I, [ pdr] = 1}, (A32) 

P:C J 



with 



S^S-f3E 



f-a- 



a2l- 



f-a- 



(A33) 



We proceed as in [42|. We first maximize S at fixed H, /, 
normalization and 'a = J prjdj]. This yields an optimal 
density pi(r,77) given by Eq. (|A30[) where (T2 is given 
by Eq. (|A3T1) -b. Then, we find that the maximization 
problem (|A32p is equivalent to 



ma_x{5[CT,C] /}, 



(A34) 



with the generalized entropy <S'[cr,^] = iS[/Cii]. An explicit 
calculation leads to 



S[a,^]=-PE^ 



tC.Q. 

02/2 . 



(A35) 



This is a sort of "mixed" case between the one studied in 
the main part of the paper (leading to the minimization 
of E'^-^-) and the one discussed at the beginning of this 
Appendix (leading to the minimization of l2^')- The 
critical points of (|A34|) return Eqs. (|A29p and (|A31|) - 
a. Furthermore, a solution of (|A32[) or (|A34[) is always 
a solution of (|A28p but the reciprocal is wrong in case 
of ensemble inequivalence. We have (|A32p <;:^ (|A34p 
(IM81) . 
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6. Conservation of E^o-, H = Hi, I = h, J?/ and 



rf-9' 



Wc consider the maximization problem 



inax{S[p]\Ef-^-, H, /, iZ/'^', [ pd?] = 1}, (A36) 
p,? J 

where S[p] is the mixing entropy (UHl) and = 
J ^prj" dydzdrj. The variations over ^ give 

/3^/' + Ai^+^2^= (A37) 
and the variations over p give the Gaussian distribution 

p= Lp-°'2v\-f^2lv\-'^ eHf^'i+o^)v ^ (A38) 



We have 



CT2 



, , (A39) 



(A40) 



2(a2 + A^2e + ^)' 



Appendix B: Detailed proof of inequality (|39p 

We consider the minimization problem 



min{ E[^,a]\H, I}, 



(Bl) 



with 



E^^J^^ dydz + \Jj dydz, (B2) 



H = S^a dydz, 



I = a dydz. 



(B3) 



(B4) 



We shall look for (local) minima of energy at fixed helicity 
and angular momentum. We proceed as in [2^, [s^] ■ The 
variations of these functionals up to second order are 



AE 



rJjS^ dydz 



SS_Sip dydz 



aSa 



■ dydz - 



{Say 



dydz, 



(B5) 



AiJ = y ^Sadydz + J a6^dydz + J S^Sadydz, (B6) 



AI = S(T dydz. 



(B7) 

The critical points satisfy the variational principle for the 
first variations 

SE + pSH + aSI = 0. (B8) 
Taking the variations over ^ and a, we obtain 

ij + pa = 0, (B9) 



2y 



(BIO) 



A minimum of energy corresponds to AE > 0. Inserting 
Eqs. dm and ((BT0| in Eq. dBS]), we find that 

AE = —p J (7(5^ dydz ^ 2 J ^^^'^ dydz 
—p [ ^Sa dydz — a [ Sa dydz ^ — [ - — ^—dydz. 

J ' J ^ J y 

(Bll) 

Then, using Eqs. dHJ and (HZ]) with AH = A/ = 0, we 
obtain 

AE — — I (5^(5 -0 dydz + I dydz + p 6 £,6 a dydz. 

2 J J J 

(B12) 

Therefore, a critical point of (jBl[) is a (local) minimum 
of energy at fixed helicity and angular momentum iff 

— J S(,Sil! dydz + J — dydz + p J S£,S<t dydz > 0, 

(B13) 

for all perturbations Sa and (5^ that conserve helicity and 
angular momentum at first order. This amounts to hav- 
ing S'^{E + pH + al) > for all the perturbations that 
conserve H and / at first order. 



Appendix C: Equivalence between (|115p and (|143p 

In Sec. I VII we have shown the equivalence of (|105p 
and (|135p for global maximization. In this Appendix, we 
show the equivalence of (jlOSp and (|135p for local max- 
imization, i.e. p(r,77) is a (local) maximum of S[p] at 
fixed E^'^', H, I and normalization if, and only if, the 
corresponding coarse-grained distribution of angular mo- 
mentum o'(r) is a (local) minimum of £''^■9- [a, ^] at fixed 
H and /. To that purpose, wc show the equivalence be- 
tween the stability criteria (|115p and (I143P. We use a 
general method similar to the one used in j22l. |42| . l49l [soj 
in related problems. 

We shall determine the optimal perturbation 5p^,{r, 77) 
that maximizes 5''J[5p] given by Eq. (|115p with the 
constraints Sa — j Sprj drj, SE-^'^' — J ipS£, dydz + 
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J Sp^ drjdydz = and J Spdr] = 0, where Sa{r) and 

S^(r) arc prescribed (they arc only ascribed to conserve 
H and / at first order) . Since the specification of Sa and 
(5^ (hence Stp) determine the second and third integrals 
in Eq. (|115p , we can write the variational problem in the 
form 



dydzdrj 



A(r)(5 / dprjdr] ) dydz 



—fl5 / Sp — dydzdrj 
4y 



C(r)<5 



where A(r), ji and C(r) are Lagrange multipliers. This 
gives 



5p*(r,7?) = -p(r,7?) 



A|^ + A(r),7 + C(r) 



(C2) 



and it is a global maximum of 6'^J[5p\ with the previ- 
ous constraints since 5^ {6^ J) = — j ^^^40^ dydzdr] < 
(the constraints are linear in dp so their second varia- 
tions vanish). The Lagrange multipliers are determined 
from the above-mentioned constraints. The constraints 
J Spdr] = and Sa = J Sprj drj lead to 



C(r) + A(r)^(r) + p 



,a2(r) 



4y 



= 0, 



(C3) 



C(r)a(r) + A(r)a2(r)4-A^ 



= -S-a{r). (C4) 



Now, the state p{r, rj) corresponds to the gaussian dis- 
tribution (|112p . Therefore, we have the well-known re- 
lations cr^(r) = CT^(r) -|- (72 and a^{r) = CT'^(r) + 3CT(r)fT2 
where (T2 = ^yj fi- Substituting these relations in Eqs. 
(jC3|) and (|C4p . and solving for A(r) and C(r), we obtain 



A(r) = -|.^(r)-A^(r), 
C(r) = A^(r)5^(r) -f- ^^'(r) - £ 



2y 



42/ 



2/3 



(C5) 



(C6) 



Therefore, the optimal perturbation (jC2[) can be rewrit- 
ten 



-|;M'7-'^) + m{:^(^--)^-^ 



(C7) 



The Lagrange multiplier p, is determined by substi- 
tuting this expression in the constraint J dydz + 

J 5/0^ dydzdr] = 0. Using the well-known identity 

cr*(r) = 5^(r) + 6a2CT^(r) + 3cr| valid for a gaussian dis- 
tribution, we obtain after some simplifications 



/i = 2/32 



ip6£,dydz+ j —Sadydz 



(C8) 



Spdr] ) dyd^; = 0, 

(CI) 



Therefore, the optimal perturbation (jC2[) is given by Eq. 
(|C7p with Eq. (|C8p . Since this perturbation maximizes 
(5^J[(5p] with the above-mentioned constraints, we have 
S^J[Sp] <S^J[Sp^]. Explicating using Eqs. ([C7)) 

and (jCSp . we obtain after simple calculations 



(5V[M < Si:SCdydz + J ^^^^ dydz 

-p j dydz — 0^ (^j ipS£,dydz + J dydz 



(C9) 



The r.h.s. returns the functional appearing in Eq. (|132p . 
We have already explained in Sec. IVI BI that for the class 
of perturbations that we consider {SH = SI = 0) the last 
term in parenthesis vanishes. Therefore, the foregoing 
inequality can be rewritten 



S'J[Sp]<^-p 



StpS^ dydz 



[Say 



2y 

p I SS^Sadydz, 



dydz 



(CIO) 



where the r.h.s. is precisely the functional appearing in 
Eq. (|143p . Furthermore, there is equality in Eq. (|C10p 
iff Sp = (5p*. This proves that the stability criteria (|115p 
and (|143p are equivalent. Indeed: (i) if inequality (|143p 
is fulfilled for all perturbations SW and (5^ that conserve 
helicity and angular momentum at first order, then ac- 
cording to Eq. (jClOp . we know that inequality (|115p is 
fulfilled for all perturbations Sp and that conserve he- 
licity, angular momentum, fine-grained energy and nor- 
malization at first order; (ii) if there exists a perturbation 
(5ct* that violates inequality (|143p . then the perturbation 
(5/9* given by Eq. (|C7p with Eq. (jCSp and Sa = Sa^, vio- 
lates (|115p . In conclusion, the stability criteria (jllSp and 
P43p are equivalent. 



Appendix D: Relaxation equations 

1. Relaxation equations associated with the 
maximization problem (|105p 

Like in classical statistical physics, it may be interest- 
ing to derive relaxation equations towards the equilib- 
rium states so as to be able to describe dynamical, non- 
stationary, regimes. On a practical point of view, these 
relaxation equations can also provide a useful numerical 
algorithm to solve the maximization problem (|105p and 
be sure that we select entropy maxima (not minima or 
saddle points). We follow the methodology described in 
[25! . We introduce a current of probability in the space of 
angular momentum fluctuations 77 and construct a set of 
relaxation equations that increase S[p] while conserving 
E^'^', I and H using a Maximum Entropy Production 
Principle (this can be viewed as the variational formula- 
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tion of Onsager's linear thermodynamics). The dynami- 
cal equations that we consider can be written as 



9C 



d 



+ V 



dp 



u- Vp 




(Dl) 



(D2) 



where X and J are two unknown quantities to be chosen 
so as to increase S[p] while conserving E^-s-, and/. In 
the second equation, the local normalization J pdrj = 1 is 
satisfied provided that J — ^ as 77 — >■ ±00. Multiplying 
Eq. (|D2[) by 77 and integrating on all the levels, we get 



9^ ^— I T 7 

— — h u • Vcr = / Jdri ■ 
ot 



Y. 



(D3) 



Next, multiplying Eq. (jD2[) by 77^ and integrating on all 
the levels, we obtain 



~dt 



u • Vcr2 = 2 / J'qdr] 



u • Vcr2 = 2 / J{ri - a)dr] 



From Eqs. ((D3)) and (|D4l) . we find that 

5(72 

The time variations of S[p] are given by 

S = - [ -TT dydzd-q, 
J pdrj 

while those of the invariants are given by 

77 



J — dydzdrj, 

2y 



= = y XtPdydz ^ 
H = = J Xadydz + J JS^dydzdrj 
1 = 0= / Jdydzdij. 



(D4) 

(D5) 

(D6) 

(D7) 
(D8) 
(D9) 



Following the Maximum Entropy Production Principle, 
we maximize S with E^-^- = H = I = and the addi- 
tional constraints 

^<Q(r,i), J ^drj<C{r,t). (DIO) 

The variational principle can be written in the form 

6S - P{t)SE^-3- - p{t)5H - a{t)6i 
1 



J i-dvj dydz 



1 



^ ] dydz = 0, (Dll) 



where /3(i), p{t), a{t), D{r,t) and time de- 

pendent Lagrange multipliers associated with the con- 
straints. This leads to the following optimal quantities 



-D 



X 



(D12) 
(D13) 



Therefore, the relaxation equation for the distribution of 
angular momentum is 



dp „ 



d_ 

drj 



D 



dp 



Integrating Eq. (jD12[) on rj we get 



Y 



D 



V 2y 



Pirn + ait) 



{D14) 



(D15) 



Inserting expressions (|D13P and (jD15[) into Eqs. (jDl[) and 
(|D3p leads to the following relaxation equations for the 
mean fiow 



9a 



u.va = -i^(^^ + /i(<)e + a(i) 



(D17) 



A relaxation equation can also be written for the centered 
variance (72. Using Eqs. (jDSl) and (|DT2|) . we obtain 



, „ 
— - + u • Vcr2 



2D 1 



2?/ 



(D18) 



Equations (|D17P and (jPlSP can be used to evaluate the 
evolution of cr^ = w'^ + <J2- The Lagrange multipliers 
evolve in time so as to satisfy the constraints. Substitut- 
ing Eqs. (|DT21) and (|DT31) in Eqs. ^7}, ([DS]) and (|D9|) . 
we obtain the algebraic equations 




it) = (D 



D 



'2y 



Pit) 

(D19) 



D 



^))m + {{x^') + {Df))pit) 

+ (D?) ait) = 0, (D20) 



) m + (D^) Pit) + ait){D) = 0. (D21) 



The coefficients D and x, which can depend on y and z, 
are not determined by the MEPP. They can be chosen 
so as to forbid divergency of the first term in the r.h.s. 
of equation (jD19|) . 

Substituting dp/ drj taken from Eq. (|D12[) in Eq. (jD6|) 
and using the constraints (jD7P - (jD9| . we easily obtain 



Dp 



dydzdrj 



X^ 



■ dydz, 



(D22) 
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so that S* > provided that D and x ^''^^ both posi- 
tive. On the other hand S = OiSJ = X = leading to 
the conditions of equihbrium and (|112|) . From Lya- 

punov's direct method, we conclude that these relaxation 
equations tend to a maximum of entropy at fixed mi- 
croscopic energy, helicity and angular momentum. Note 
that during the relaxation process, the distribution of an- 
gular momentum is not Gaussian but changes with time 
according to Ec^. (jD14|) . The distribution is Gaussian 
only at equilibrium. Therefore, these relaxation equa- 
tions describe not only the evolution of the mean flow 
but also the evolution of the distribution of fluctuations. 
We stress, however, that these equations are purely phe- 
nomenological and that there is no compelling reason 
why they should give an accurate description of the real 
dynamics. However, they can be used as a numerical al- 
gorithm to compute the equilibrium state corresponding 
to (|105p . Indeed, these equations can only relax towards 
an entropy maximum at fixed microscopic energy, helic- 
ity and angular momentum, not towards a minimum or 
a saddle point that are linearly unstable with respect to 
these equations. 

Remark: In fact, we will find in [2lj that there is no 
entropy maximum, just saddle points. In that case, the 
dynamical equations lead to a "collapse" at smaller and 
smaller scales, similar to the Richardson energy cascade 
in 3D turbulence. However, we will also observe that 
the system can remain blocked in a large-scale coherent 
structure (like in 2D turbulence). In the present 2.5D 
situation, this is an unstable state (saddle point of en- 
tropy) , but it can persist for a long time if the dynamics 
does not spontaneously develop the "dangerous" pertur- 
bations that destabilize it. This is because a saddle point 
is unstable only for some perturbations but not for any 
perturbation. 



where f3{t) is determined by the constraint on the micro- 
scopic energy leading to 



^ = 2Ef '- - y dydz - y ^ dydz. (D26) 
On the other hand, the time variations of H and / are 
H = = J X a dydz + j Yldydz, (D27) 
i = 0= Ydydz. (D28) 



Following the Maximum Entropy Production Princi- 
ple, we maximize S with I = H ~ (the conservation 
of the microscopic energy has been taken into account in 
Eq. (|D26|) ) and the additional constraints 

^<Q(r,i), ^<C.(r,t). (D29) 
The variational principle can be written in the form 

dS - ii{t)5H - a{t)5i 

1 f X'^\ 

-5 I -— I dydz 



X(r,t) V 2 
6 —— dydz = 0, 



D{v,t) V 2 
and it leads to the following quantities 

x = -x m)^ + Km 



Y ^-D 



V 22/ 



(D30) 

(D31) 
(D32) 



Inserting expressions (|D3ip and (jD32[) into Eqs. (jD23[) 
and (|D24[) . we obtain the relaxation equations 



Relaxation equations associated with the 
mciximization problem (|124p 



We shall now introduce a set of relaxation equations 
associated with the maximization problem (|124p . We 
write the dynamical equations as [53| : 



dt 



u - = 



d_ 

dz 



V 



X. 



(D23) 
(D24) 



^+u-Va^-Di^^+f^m + ait)y (D34) 

The Lagrange multipliers evolve so as to satisfy the con- 
strain ts. Substi tuting Eqs. (|D3ip and (|D32p in Eqs. 
(|D27)) and (|D28l) . and recalling Eq. (|D26| . we obtain the 
algebraic equations 

^^2£'* -(?«)- (I), (D35) 



where X and Y are two unknown quantities, to be chosen 
so as to increase S'[^,a'] while conserving E^s-, H and / 
given by Eqs. (fT25)) . ((^61) and (|T27)) . The time variations 
of S are 



+ (ZJj) a{t) = 0, (D36) 



-i^(t) ^2 J i)Xdydz + j ^Y dydz ] , (D25) 



D- ) m + {DO Kt) + ait){D) = 0. (D37) 
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Substituting ijj and a/y taken from Eqs. ()D31I) and (|D32[) 
in Eq. (|D25P and using the constraints (|D27p and (|D28p . 
we easily obtain 



S 



X 



■ dydz 



D 



dydz, 



(D38) 



so that S' > provided that D and x are both positive. 
On the other hand 5' = 0ifFX = r = leading to the 
conditions of equilibrium given by Eqs. (|130|) and (|13ip . 
From Lyapunov's direct method, we conclude that these 
relaxation equations tend to a maximum of entropy at 
fixed microscopic energy, helicity and angular momen- 
tum. 

The relaxation equations (jD33P and (jD34[) arc similar 
to Eqs. (jD16p and (jD17[) but the constraints determining 
the evolution of the Lagrange multipliers are different. 
More precisely, Eqs. (jD36|) and ()D37[) coincide with Eqs. 
(|D20]) and ((MT|) but Eq. (|DT9|) has been replaced by 
Eq. (|D35p . Indeed, in the present approach, the distri- 
bution of angular momentum is always Gaussian during 
the dynamical evolution. It is given by Eq. (|117p at any 
time, i.e. 



(D39) 



\ Airy J 

By contrast, in Sec. IVI Al the distribution of angular 
momentum changes with time. Therefore, the dynamical 
evolution is different. However, in the two approaches, 
the equilibrium state is the same, i.e. it solves the maxi- 
mization problem (|105p . This is sufficient if we use these 
relaxation equations as numerical algorithms to compute 
the maximum entropy state. 

Remark: Using Eqs. (|D23l) . (|D24l) and (|D25|) . it is 
easy to show that S = -l3{t)E''-s- so that E^'-s- < since 
/3{t) > 0. Therefore, the macroscopic energy monotoni- 
cally decreases through the relaxation equations. This is 
to be expected since the maximization problem (|124p is 
equivalent to the minimization of the macroscopic energy 
at fixed helicity and angular momentum (see Sec. IVI CI) . 



3. Relaxation equations associated with the 
minimization problem (|135p 

We shall introduce a set of relaxation equations asso- 
ciated with the minimization problem (jl35p . We write 
the dynamical equations as 



X, 



(D40) 
(D41) 



where X and Y are two unknown quantities to be chosen 
so as to decrease E'^-^- while conserving H and /. The 
time variations of iJ'^-s- are given by 



ipX dydz - 



2y 



Y dydz. 



(D42) 



On the other hand, the time variations of H and / are 
H^O^ J XWdydz + J Y^dydz, (D43) 
i = 0= I Ydydz. (D44) 



Following the Maximum Entropy Production Princi- 
ple, we maximize the dissipation with I = H = 
and the additional constraints 

4^<Q(r,0, ^<a(r,i). (D45) 
The variational principle can be written in the form 
SE^-s- + n{t)SH + a{t)6i 

1 ^ x^\ 

-S —— ) dydz 



X(r,t) V 2 
-6 -— dydz ~ 0, 



(D46) 



and wc obtain the following quantities 

X = -X + ix[m , (D47) 

(D48) 



Y = -D[-+^,{t)i + a{t) 



Subst ituting Eqs. (jPTT]) and (jD48|) into Eq. (|D401) and 
(jD4ip leads to the following relaxation equations 

| + u.Ve=|(J)^x(^ + .W^), (D49) 

^ + u.Va = -i^(^^-f^(i)? + a(t)) . (D50) 

The Lagrange multipliers /i(t) and a{t) evolve so as to 
satisfy the constraints (|D43p and (|D44p . Substituting 
Eqs. (|D47| and (fPiSl) in Eqs. (|D43l) and (jD44| . we 
obtain the algebraic equations 

+ {xi>a) + (d^^ = 0, (D51) 



{DC)n{t) + o^it){D) + {D—)^{). (D52) 



Substituting ip and a/y taken from Eqs. (|D47j) and (|D48j) 
in Eq. (|D42|) and using the constraints (|D43|) and (|D44p . 
we easily obtain 

E'-a- = - / dydz ~ / — dydz, (D53) 

J X J D 

so that E'^-3- < provided that D and x are both posi- 
tive. On the other hand, E^'J- ==OiffX = y = leading 
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to the conditions of equilibrium given by Eqs. (|14ip and 
p42p . By Lyapunov's direct method, we conclude that 
these relaxation equations tend to a minimum of macro- 
scopic energy iJ'^f - at fixed helicity and angular momen- 
tum. Therefore, the relaxation equations (|D49ID50[) can 
be used as a numerical algorithm to solve the minimiza- 
tion problem (|135p . 

Remark: since these relaxation equations solve Eq. 
((32)) . they can also be used as a numerical algorithm to 
construct nonlinearly dynamically stable stationary solu- 
tions of the axisymmetric Euler equations corresponding 
to Beltrami states (see Sees. IIIII and ITV ]) independently 
of the statistical mechanics interpretation. 



Following the Maximum Entropy Production Princi- 
ple, we maximize the dissipation with H = and 
the additional constraints 



(E5) 



The variational principle can be written in the form 



(E6) 



Appendix E: Another type of relaxation equations 

In the main part of the paper, we have not taken into 
account the conservation of circulation ^ = J S, dydz be- 
cause there is no critical point of energy at fixed helicity, 
angular momentum and circulation (see [2l|). Neverthe- 
less, at the level of the relaxation equations, it is possible 
to take this constraint into account. We shall introduce 
a set of relaxation equations that minimize the energy 
iJ'^'S' at fixed helicity H ^ angular momentum / and cir- 
culation r. Since there is no energy minimum (not even 
a critical point of energy), these equations should have 
a non-trivial behavior. To derive these equations, one 
possibility is to write them in the form (|D40p - (jD41[) and 
introduce Lagrange multipliers for each constraint. An- 
other possibility is to write them in the form 



d 



dt 



dz \Ay 



da 
'dt 



u • Vct 



-V- J, 



(El) 
(E2) 



where and Jo- are two unknown currents to be chosen 
so as to decrease E'^-^- while conserving H. The form 
(|Eip - (|E2p guarantees the conservation of circulation and 
angular momentum. The time variations of i^^f- are 
given by 



E^-s- = y J? • yi'dydz + y • V i^— j dydz. (E3) 
On the other hand, the time variations of H are 



and we obtain the optimal currents 
J. = -D. 



( W + /i(t)VCT) 



V( - ) +M(tm 



(E7) 
(E8) 



Substituting Eqs. dET]) and dES]) into Eq. (|ET]) and (|E2 
leads to the following relaxation equations 



da 
'dt 



u • Va 



V ( - ) + /i(t)Ve 



(E9) 
(ElO) 



The Lagrange multiplier /i(t) evolves so as to satisfy the 
constraint (pij) . Substituting Eqs. ([E7| and (pS]) in Eq. 
(p4)) . we obtain 



^l{t) 



f D^yij ■^adydz + J D„\/ (^) -S/^dyd. 



J D^{VaY dydz + J D„{V£,Y dydz 



(Ell) 



Substituting Vip and V{a/y) taken from Eqs. (jE7p and 
(jESP in Eq. ()E3P and using the constraint (|E4p . we easily 
obtain 



E^-s- - - / ^ dyd^ " / ^ '^y'^^' ^^^^^ 



H = 0= / • V a dydz + V^dydz. (E4) 



so that E'^-^- < provided that and Do- are both 
positive. 
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